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Abstract 


A basic problem, which is of considerable interest in 
geoastrophysical applications of magnetogasdynamics , is the 
nonlinear interaction of a fast shock (S^) with a tangential 
discontinuity (T) . The problem is treated for an arbitrary 
interacting with an arbitrary T under the assumption that in 
the frame of reference in which and T are at rest the flow 
is super-fast on both sides of T and that a steady flow develops. 
As a result of the nonlinear analysis a flow pattern is obtained 
consisting of the incident discontinuities S^l and T2 and a 
transmitted fast shock S^3, the modified tangential discontinuity 
T4 and a reflected fast shock S^5 or fast rarefaction wave R^5. 
The results can be discussed in terms of seven significant 
similarity parameters. The uniqueness of the solution is in- 
vestigated as well as the region of dimensionless parameter 
space in which a solution of the type described above is 
possible. Varying the similarity parameters the solutions of 
type S f 1,T2,S f 3,T4,S f 5 are bounded in parameter space by the 
disappearance of S^5j by the maximum total pressure ratio across 
S f 5 at which the flow direction on both sides of T4 can just be 
matched and by the limiting cases of S f 3 or S f 5 being switch-on 
shocks. The following analytical and numerical results concern- 
ing the dependence of the interaction on the type of tangential 
discontinuity T2 are obtained: A T2 with a velocity shear only 
in the direction perpendicular to the plane XY spanned by the 
normals of S^I and T2 leaves the incident shock unchanged and 
no reflected wave S^5 or P^5 develops. A velocity shear across 
T2 in the plane XY can be compensated by a density change such 
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that the incident shock is not deflected or changed at all 
and that no S^5 or R^5 develops although appreciable changes 
of density and velocity may occur across T2. In addition 
special cases like changes in magnetic field direction only, 
changes in density or velocity shear only etc. are discussed 
in some detail. 
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I . Introduction 


In recent years plasma and magnetic field experiments on 
satellites and space probes have proved the existence of all 
proposed types of magnetogasdynamic discontinuities in the in- 
terplanetary plasma with the exception of the contact discon- 
tinuity with nonvanishing normal magnetic field. A comprehen- 
sive classification of the possible discontinuities taking into 
account the evolutionaritv conditions is given in the book by 
Jeffney and Taniuti . The pertinent space observations are 
reviewed in a recent article by Burlaga^. For many problems of 
the physics of discontinuities in the collisionless interplane- 
tary plasma, magnetogasdynamics or even gasdynamics can give a 
reasonable answer, although it is realized that the properties 

of the discontinuities are modified by pressure aniso- 
3 4 5 

tropies ’ ' , the possible change in heat flux across the dis- 
continuity 6 and a background wave spectrum 6 ' 7 . 

The observations show more specifically that in addition 
to continuous fluctuations of the physical quantities the 
spatial and time structure of the solar wind is characterized 
by the frequent occurrence of tangential discontinuities and 
in somewhat smaller numbers of rotational discontinuities 8 . On 
the other hand, the interplanetary plasma is the propagation 
medium for fast magnetogasdynamics shocks which produce major 
disturbances in the space, environment of the earth and there- 
fore are of considerable interest. 

A basic physical process during the propagation of inter- 
planetary shocks is their interaction with the large number of 
tangential discontinuities . Other applications of this basic 
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process are the interaction of tangential discontinuities with 
the earth’s bow shock and the interaction of interplanetary 
shocks with the magnetopause of the earth after modification 
of the shock by its collision with the earth's bow shock. 

It is the purpose of this paper to investigate the non- 
linear interaction of a fast magnetogasdynamic shock wave with 
a tangential discontinuity in the nonrelativistic case. Both 
discontinuities are considered to be arbitrary in their orien- 
tation to each other and in all other physical parameters with 
the only restriction being that in the rest frame of both dis- 
continuities a steady flow pattern is possible. Other possible 
solutions include a break-up of the shock, a propagation of 
part of a disturbance upstream along the initial discontinuities 
etc. The case of a steady flow is most important because the 
flow pattern depends on the local physical properties in the 
vicinity of the line of intersection of both discontinuities 
only. 

In the discussion of the numerical results we shall empha- 
size the case of only fast shocks in the resulting flow. Two 
limiting cases of this more general problem have been considered 
so far. Jaggi and Wolf 9 have investigated the problem of a weak 
tangential discontinuity interacting with a fast shock. The case 

with vanishing magnetic field is well known in the gasdynamic 
1 0 

literature . The special case, where both discontinuity planes 

are parallel, is included in the discussion of the general mag- 

1 1 

netogasdynamic Riemann problem by Gogosov 

In Section II we derive the. resulting flow pattern from 
general arguments. In Section III the equations for the solution 
of the problem are formulated and seven dimensionless parameters 
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are derived, which characterize the problem and lead to some 
interesting similarity laws. Section IV contains a graphical 
representation of the mathematical problem and a discussion 
of the uniqueness of the results. In addition, those regions 
in parameter space are identified in which solutions of the 
type considered are possible. 

In Section V some numerical results for the interaction 
are presented. 

Except for the requirement that each of the participating 
discontinuities has to be evolutionary, we shall not consider 
the stability of the solutions. 
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II. Derivation of Basic Flow Pattern 


The problem under consideration is the interaction of a 
plane fast magnetogasdynamic shock and a plane tangential dis- 
continuity, which according to the nomenclature of Jeffrey and. 

12 

Taniuti is a contact discontinuity with a vanishing magnetic 
field component in the normal direction. Both discontinuities 
may have an arbitrary orientation to each other. If the planes 
are not parallel, it is always possible to find the line of 
intersection of both discontinuities, which we shall call the 
interaction line. We can then always find a frame of reference 
in which the fast shock and the tangential discontinuity are at 
rest. In this frame of reference we define a right-handed 
Cartesian coordinate system X,Y,Z with the Z-axis parallel to 
the interaction line and the X-axis in the direction of 
propagation. This is shown in Fig. 1 , which will be used to 
develop the basic magnetogasdynamic flow pattern throughout 
this section. The incident fast shock is shown as a solid line 
called S^I and the initial tangential discontinuity T2 is 
shown as a dashed line. A straight line in the plane of 
Figure 1 evidently corresponds to a plane perpendicular to the 
XY plane, since the problem is two-dimensional. In defining 
our frame of reference we have not yet made use of the degree 
of freedom in the Z-direction. For definiteness we assume the 
velocity component v . in sector 1 between S-1 and T2 to be 

u f * r 

zero without loss of generality. 

The physical state in sector 1 characterized by the ve- 
locity vector v^ , the magnetic field vector , the pressure p^ 
and density is given as well as the physical quantities , 
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^ 2 '’ P 2 " p 2 on t ^ ie ot ^ er s ^ e °f T2. It is clear that it is 
impossible to find a steady-state solution for all choices of 
initial parameters. The parameter range for which this is 
possible will also be investigated in this study. 

It is illuminating to consider the cases in which one of 

the initial discontinuities 5^1 and T2 is weak, where the 

weakness is characterized by a parameter e << 1 . The case of 
9 

T2 being weak leads to a sector pattern, in which the shock 
S^I is changed by an amount of 0(e) resulting in a transmitted 
fast shock S^3. T2 is modified into the tangential discontinuity 
T4 with a strength of 0(e) . A reflected weak fast shock S^5 or 
rarefaction wave R f 5 is generated with the strength of 0(e). 

The direction of S^5 or R^5 is given by the fast Mach line 
going out from the interaction line and making an acute angle 
with the Y-axis. Here we have used the following convenient 
notation. The homogeneous sectors between the various dis- 
turbances are counted counterclockwise giving the subscripts 
for the physical quantities in these sectors. The disturbances 
are counted in the same sense and the characters R, S, T give 
the type of discontinuity. Physical quantities defining the 
disturbances themselves have subscripts according to this 
counting. The subscripts f and s define fast or slow shocks (S) 
or rarefaction waves (R) in the magnetogasdvnamic sense. 

The results for T2 being weak of 0(e) suggest by the 
continuous dependence of the solutions on the initial conditions 
that for a range of tangential discontinuities T2, which are 
finite in amplitude, the flow pattern will be S f 1 , T2, S f 3, T4, 
R f 5 or S^5. For the appropriate range of velocities there 
could also be real characteristics or Mach lines of the slow 
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mode in sector 2 and therefore solutions involving slow shocks. 
However, we shall restrict ourselves to the case of fast trans- 
mitted waves only. 

These considerations could be used to establish the basic 
magnetogasdynamic flow pattern in the frame of reference chosen 
by their extension to the general case of S^I and T2 being of 
finite strength. A more satisfactory and instructive way is to 
derive this basic flow pattern from basic principles and the 
known properties of magnetogasdynamic discontinuities and rare- 
faction waves. This is especially important for the question of 
the uniqueness of the solutions. 

We derive the basic type of solution by the following five 
arguments . 

1 . Since the initial conditions given by the angle a 2 (see 

Fig. 1) and the physical quantities v^ , , p^ r p^ and v 2 , B 2 , 

p 2' p 2 do not contain any characteristic length and since the 

magnetogasdynamic equations contain no characteristic length, 

we immediately get the result that the solution to our problem 

depends on the angle cf only and not on r, where r, cf> and Z are 

£ 

cylindrical coordinates around the Z-axis, i.e., we have — = 0 

3 9r 

and — = 0. The solution must then be composed of a sector 

3Z 

pattern of discontinuities and rarefaction waves centered at 
the origin, and homogeneous regions in between. 

2. The vector v^ and must be parallel to the plane of 
T2 by the definition of a tangential discontinuity, i.e., the 
projections of v 1 and B. on the X, Y-plane called v . and 

it Pf i 

Bp i , respectively, must be parallel or antiparallel to each 
other. Because of E = -v X B we must therefore have = 0 in 
sectors 1 and 2 and because of Faraday's law = 0 everywhere. 
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Therefore v p and B p must be parallel or antiparallel every- 
where. This immediately excludes a contact discontinuity with 
nonvanishing normal component of the magnetic field B n ^ 0 from 
further consideration, because it would imply = - (v X B) z # 0. 

3. There can be at most one tangential discontinuity in 
addition to T2 as the result of the interaction. Suppose, there 
are more than one newly created tangential discontinuities among 
them T^ and T n with n > k > 2, say. Since a tangential dis- 
continuity is characterized by vanishina normal components of 

B 2 

v and B and the continuity of the total pressure p t = p + g^- 
across the discontinuity only, the state between T k and T n is 
not completely determined by the initial conditions. There is 
for example no way to determine the density or densities between 
T k and T . As we have to require at least piecewise continuous 
dependence of the solutions on the initial conditions such a 
configuration is impossible unless it is required by additional 
initial conditions. The statement at the beginning of this 
paragraph is proven thereby. 

4. An important quantity in the classification of dis- 
continuities is the Alfven Mach number A based on the normal 

n 

components of the velocity and magnetic field with respect to 
a shock or rarefaction wave under consideration. Using the 
normal vector n, xvhich points in the direction of propagation 
of the shock or rarefaction wave, we obtain 

A = -I L v/4irp = -* — 2 L /4irp = ‘ •/ 4 it p ( 1 ) 

|B-n| iB p =n| | B p | 

since the direction cosines of v p and B p with respect to an 
arbitrary normal n are equal in magnitude as shown under item 2 . 
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That is, the Alfven Mach number A n is independent of the normal 
direction, if the state in front of the wave is specified. This 
property has interesting consequences. Before we investigate 
these consequences some inequalities for shocks and rarefaction 
waves (including strength zero) have to be recalled! 

v 

S„s upstream; A >1, — > 1 (2) 

r n — c^ — 

v 

downstream; A > 1 , — < 1 
n — ' c f — 


v 

R^; upstream and downstream; A^ >_ 1 , — 

v 

S ; upstream; A > 1 , — > 1 
s n — c — 

s 

v 

downstream; A > 1 , — < 1 
n - c s - 

v 

R ; upstream and downstream: A < 1 , — 
s r n — c 

s 


1 


(3) 

(4) 

(5) 


Here c^ and c s are the fast and slow magnetoacoustic speeds, 
respectively. For finite shocks A n = 1 is valid on the down- 
stream side for switch-on shocks only and on the upstream side 
for switch-off shocks only. 

Using equation 1 and inequalities 2, 3, 4, 5 it immediately 
turns out that a streamline starting in sector 1 can pass through 
fast waves (i.e., fast shocks or rarefaction waves) only, because 
the initial shock wave has been chosen to be of the fast 

type, i.e., A n > 1 everywhere on the streamline starting in 
sector 1 . Here we always assume that no switch-on or switch-off 
waves occur as part of the general solution. We shall return to 
this question later. A stream line starting in sector 2 will 
also pass through 'fast waves only, if 

h „ > max (1 /s~) 
n , 2 — 2 
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( 6 ) 



(7) 


where S 2 is given by 


s 2 = 5/3 


P2 

b2/4tt 
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as the square of the ratio of sound-speed over Alfven speed ' 

It has already been mentioned that we shall consider the case 
of fast transmitted waves only which requires inequality 6 to 
be valid. 

5 . We now have to introduce the important notion of 
"outgoing'' waves and "incoming" waves, where by waves we mean 
shock waves or rarefaction waves. "Outgoing" fast waves are 
defined such that a weak fast disturbance can travel along the 
wave away from the interaction line and thereby modify the wave. 
It is based on the requirement that a change in the initial 
conditions must influence the various disturbances resulting 
in the S f 1 , T2 interaction. This implies, that behind the shock 
or downstream of the shock the state must be sub-fast, super- 
Alfvenic and super-slow. This condition restricts the allowable 
angle with the streamlines of the incident plasma as schemati- 
cally shown in Figure 2a. Since the fan which characterizes a 
fast rarefaction wave is made out of fast characteristics, we 
have the situation shown in Figure 2b. We shall call "incoming" 
waves those waves which are not "outgoing". 

Since in the interaction only outgoing waves can be cre- 
ated, i.e., disturbances which are connected to the interaction 
line by causality, the number of possible shocks or rarefaction 
waves going out from the interaction line can be further reduced. 

For example, behind an outgoing fast shock or rarefaction wave 
v n - - 

A n > 1 and — < 1 for all possible normal directions of a wave 
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following the outgoing fast wave. Therefore, no further shock 
or rarefaction wave is possible on a streamline passing 
through an outgoing fast shock. 

A streamline starting in sector 2 can therefore pass 

through one outgoing fast wave S^3 or R^3 only and has to reach 

infinity in sector 3. If S^l is an outgoing shock no further 

fast wave could be crossed by a streamline until infinity is 

reached behind S^I . This is impossible, since there is only one 

unknown namely the strength of S^3 or R.^3, whereas two conditions 

have to be fulfilled when sector 3 and the state behind S f 1 are 

separated by a tangential discontinuity or seven conditions if 

all physical quantities have to be continuous from sector 3 to 

the state behind S X 1 . Here we have used v X B =0 everywhere. 

f -P -P 

Consequently S f 1 must be an incoming shock in agreement 
with the fact that S f 1 is given by the initial conditions. The 
general solution can then be given as S f 1 , T2 leading to S f 3, 

T4, R^5 or S^5. We also have five homogeneous sectors 1 to 5, 
where 5 is the sector behind S^l . With the arbitrary strength 
parameters e 3 and e 5 for S f 3 and S f sor R^5 considered as un- 
knowns, two equations have to be fulfilled 

p t,3 * e 3^ = p t,4 ^ 

and a Vf3 (e 3 ) = V 4 U 5 ) < 9 > 


where a v ^ is the angle between v^ and the Y-axis in sector j . 
A' solution involving a rarefaction wave R^3 is not possible, 
because it would require a v 3 < a v 4 • Because of p fc 5 > p t ^ , 
equation 8 and 1 = Pt 2 we ^ et p t 5 > p t 4 ' wh:i - cl1 requires 


12 



and 


a wave R f 5. Since a v ^ = a v>2 and a v , 4 > 
a , < a ,, we obtain a - < a - in contradiction to 
equation 9 . 

The condition, that S^I is an "incoming" wave, leads to a 
limitation of 


“M f ‘ “2 < “crit * 90 


( 10 ) 


where a., is the Mach angle of the fast wave i.e., the angle 

Mr 

f 

a 2 at which the fast Mach number equals one for given s^ , 

a . _ r and A ... The angle a . . is given by solving equation 
io^ i n, I 

1/2 


.2 1 S 5 + 1 . //i S 5_ 

n,S" 1 “72 l\2 77X 

c °s a iB>5 \ cos rf iB>5 


14 - 

K/ 


2 

COS a 


iB,5> 


( 11 ) 


for a 2 = a cr it' where A n 5 and s 5 are functions of , a 2 , 
a._ - and A - and where s., a.„ - and A 1 are considered to 
be given, is defined analogously to s 2 in equation 7. a ig ^ 
is the inclination of the magnetic field in sector 1 with re- 
spect to the XY-plane with a ig 1 = O for B z ^ = O and 


B y ,1 > °' 
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Ill . Similarity Laws 


We 


rewrite equations 8 and 9 using P t = P-t 2 an< ^ 


a . = a _ 
v, 1 v, 2 


p t,3 P t , 4 P t , 5 


or 


Y t,f (s 2 7 a iB,2' A n,2 7 a 3 " a v , 2 5 ~ Y t , f (s 5 7 a iB,5 7 A n,5 ,a 5 + “v,5 ) 


Y t,f (S 1 7 a iB , 1 7 A n,1 7 a 2 ) 


( 12 ) 


and 


a _ ~ a _ = a . - a c + a c - a ,, 

v,3 v,2 v,4 v,5 v,5 v,1 


or 


V s 2' a iB , 2 ' A n,2' “3 ' “v,2 )= - 5 f (s S’ a iB,5' A n,S' V a v,5> 


+ « f (Si> a iB , 1 ' A n,1' “2* (13) 

Equations 12 and 13 have to be solved for the unknowns a 3 and 
a 5 „ The angles a 3 and a 5 are the angles between the dis- 
continuities S f 3 and S f 5 and the Y-axis as shown in Fig. 1. 

R f 5 is characterized by the angle between the weak discon- 
tinuity at its downstream side and the Y-axis. The functions 
Y f and are derived from the general jump relations in the 
appendix for shocks. The appendix also contains the functional 
dependence of s 5 , “ iBf5 , A n , 5 ? - nd <* v?5 on s 1 , « 2 , o iB ^, A^. 

Inspection of equations 12 and 13 shows, that the initial 
conditions define seven dimensionless parameters s 1 , a 2 , a. B ^, 

A s~,a._ 0 and A 0 , which characterize the complete problem, 
n, 1 2' iB,2 n,z 
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We have started with 17 quantities p ^ , p^ , , P 2 , P 2 ' 

y_ 2 , B 2 or 14 quantities, if the jump relations across T2 are 
taken into account. This reduction in parameters to seven di- 
mensionless parameters is very convenient for a further dis- 
cussion and leads to what can be referred to as similarity 
laws. The problem of the uniqueness of the solutions of 
equations 12 and 13 will be treated in Section IV. 

The solution determines the dimensionless quantities 
characterizing each sector and each discontinuity or rarefaction 
wave, i.e. , o^, , a iB/ j, a v ,j' A n,j' r B,j with 3 = 3, 4, 5 and 

the fast Mach numbers for S^3 and S^5 or R f 5, where r fi ^ is 
the magnetic field magnitude in sector j measured in units of 



The parameters entering equations 1 2 and 1 3 are independent 
of the velocity component in the Z-direction. Introducing a iv ^ 
the inclination of the velocity vector in sector j the quanti- 
ties a. and a. . are determined by the choice a. . = 0°. 

To determine a. it is necessary to specify an eigth initial 
parameter ou v 2 „ a^ v 2 only influences a iv? 3 , however. 

Obviously, the solution for given physical parameters c* 2 , 
p 1 ' Pi / , B-j , p 2 , P 2 v 2 , B 2 can easily be obtained by using 

the results of the dimensionless problem. The initial conditions 
must be compatible with a tangential discontinuity connecting 
sector 1 and sector 2. 

Symmetries in the discontinuity jump relations immediately 


lead to a number of simple results . Replacing a by - a . 

Id / I < IB f \ 


leads to a 


a iB, 5 ’ a iB,4 


a iB , 4 * a iB,2 " “ a iB,2 in 


sector 2 leads to a._ _ -*■ - 

1 D f J 


a._ _ in sector 3, etc, 

^ J 
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Jaggi and Wolf consider the generation of fast magneto- 
acoustic waves (in our notation S ^5 or F.^ 5 ) by weak tangential 
discontinuities % They show results for all five linearly inde- 
pendent types of discontinuities. As our dimensionless para- 
meters show, there are only three parameters s, a.~ and A , the 
changes of which characterize the kind of interaction generating 
S^5 or R^5. Their Figures 5, 7, 8 can therefore be simply re- 
lated by comparing the corresponding changes in A n from sector 1 
to sector 2 . 

The similarity laws included in the seven or eight di- 
mensionless parameters can be used to immediately derive some 
important results. Suppose a 2 < a cr; j^ t° given together with 

P-|o P 1 » v 1 , B 1 and P 2 = P 1 ' p 2 = P 1 ' -2 = -1 and B 2 = . Since 

T2 is a zero-strength discontinuity we would expect no inter- 
action and no change in the shock flow pattern at all, i.e. 

S f 3 would be identical to S f 1 and there would be no S f 5 or R f 5. 
There would be a zero-strength discontinuity T4 with an angle 
“4 < « 2 becaus e of the bending of the stream-lines as they pass 
through S f 1. wo also have s 2 = Sl , <» iB(2 = « 1B>1 , A „ /2 = 
and “iv ,2 = a i V/ 1 * If we introduce a velocity shear in the 
Z-direction by choosing v^ 2 * 0, none of the seven charac- 
teristic parameters, and therefore also the solution just de- 
scribed, change except for v Zf3 = v z 2 + v z 5* Consequently, a 
shock S f 1 interacting with a tangential discontinuity with a 
velocity shear in the direction of the interaction line only 

will not be modified, i.e., n- , = n„ and all properties of 

L 'f f 1 

the shock S f 1 are identical to those of shock S f 3. Only the Z- 
components of the velocity vector will be different for both 
shocks. The angle oi 4 will be less than a 2 and T4 is characterised 
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by the same velocity shear as T2 and no change in any other 
physical quantity. There will be no generated or reflected 
wave Sj5 or R^5. The same nondimensional solution applies for 
T2 having a velocity shear in the XY-plane, if it is compensated 
by a density change such that A. = A « and s~ = s. as well as 
a iB 2 = a iB 1 * For the densit Y ratios we must have 

2 2 
Po 1 p_ v. 

_2 = Pf 1 _ _3 _ p ; 4 

" ^.r -.4 ' ■ 

These considerations show that a small change in the di- 
mensionless quantities s, A and ct . _ across T2 is a sufficient 

n I# 

condition for the interaction to be weak, i.e. f for the change 
S f 1 S^3 to be small and the strength of S^5 or R^5 to be weak. 
This condition is evidently much less restrictive than the 
condition of a small change in the dimensional quantities p, p, 
v, B. 
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IV. Uniqueness of Solution and Boundaries in Parameter Space 


The uniqueness of the solutions and can be discussed 
with the help of the graphical representation in Figure 3 in 
the a v , P t /P t -j plane. We consider the seven initial para- 
meters s 1# a 2 , a ±Bf1 , A n ^, s 2 , a iB f 2 ' A n ,2 t0 be given " The 
point Se 1,2 then gives the state in sectors 1 and 2. The curve 


labelled S-1 is the locus of all points (a - , p _/p 1 ) 

obtained by keepinq s. , a.., . and A w in sector 1 constant 

I XJD g I n g I 

and varying a 2 between (s 1 , a iB ^ 1# A n 1 ) and a crit (s 1 , 

a iB,V ^n 1 ^ • BG Y° nd the point , corresponding to 
“2 = a crit' tbe s h°ck wave S^I would be outgoing. The initial 
value of a 2 chosen in the initial parameter set leads to point 
Se5, which is also the upstream state for S f 5 or R f 5. The curve 
labelled S f 5 is the locus of all points (a . , p ./p ) 

obtained by varying « 5 such that a shock S f 5 is given with the 

upstream state s 5 , a iB 5 , A fi It contains outgoing and in- 
coming solutions. The outgoing and therefore allowed part of 
the curve S^5 starts at Se5, passes through Se4, T and Se'4 and 
ends shortly thereafter. The curve S^5 is connected continu- 
ously up to second order derivatives with the rarefaction wave 


branches R f 5, which are not shown completely. Similarly, the 
curve S^3 is the locus of the S^3 shocks obtained by varying 
a 3 - It depends on s 2 , a iB , 2 and A n 2 ° The intersect i°ns of 
curve S^3 with curve Sj5,R^5 give the solutions to equations 8 
and 9 or 1 2 and 13. Every point on these curves is uniquely 
given by an angle a 3 or a 5 , respectively. The case shown 
permits two solutions given by the points Se4 and Se'4, in 
which S f 3 and S f 5 are both "outgoing". However, the curves S f 3 
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and S^5 could also include switch-on shock waves for s 2 < 1 
and/or s,. < 1 , respectively. For example, the curve 3^5 then' 
ends in a switch-on shock as a limiting case, before it reaches 
the line = a v 5 with infinite slope. In such cases there is 

the possibility of one solution only. If S^3 misses the allowed 
part of R^5, S^5, no solution is possible. One invariable feature 
of these curves is that the curve labelled S f 5, R f 5 is symmetric 
to the horizontal line a v 4 = a v 5 / if for a moment we do not 
distinguish between "incoming" and "outgoing" waves. In addition, 
on each of the symmetric parts the second derivative does not 
change sign. 

If more than one solution is possible, we pick the ad- 
missible one by the fo3.1owing postulates The admissible solution 
Se4 must continuously tend to the point Se5 without S f 3 ever 
being tangent to S^5 or R^5 during the approach, when s 2 , 2 

and A n ^ 2 tend to , a iB ; 1 ' A n, 1 ' res P ectivel Y / in a continuous 
manner. Although this device works remarkably well in gas- 
dynamics, a more physical selection of the correct solution 
would be desirable. It could possibly be provided by a stability 
analysis of the type carried through for a plane shock by Gardner 
and Kruskal 14 . The solution of this problem is outside the scope 
of this work. It would be interesting to investigate experimen— 
tally the question whether under certain conditions a flipping 
from solution Se4 to SeM is possible as a result of an in- 
stability of solution Se4 . Solutions of the type SeM could 
also occur as part of a more complicated interaction pattern 
outside the range of the initial parameters making a steady 
flow pattern possible. The numerical results in Section V will 
only contain the admissible solutions Se4. Finally we note a 
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necessary condition for Se4 to be admissible. It may be ad- 
missible only, if it lies on the branch from Se5 to the point T 
in Fig. 3, which is obtained by drawing the tangent from Se1,2 
at curve S^5. 

An important problem is to find the region in the seven- 
dimensional space of our initial parameters, in which the 
steady-state solution is possible. If the orientation of both 
discontinuities and the physical parameters in front of S^l 
are given, s^ , a 2 , A n 1 are determined. From Figure 3 

and the previous discussion we see that three limiting cases 
exist for the basic steady flow pattern S f 1 , T2, S f 3, T4, S f 5. 
Firstly the shock S^5 can tend to zero strength, i.e., the 
fast Mach number M f , 5 = 1. Secondly, the curves S f 5 and S f 3 
can become tangent. After a further change of the parameters 

s 2' a iB,2 ? A n,2 r ^^ht direction, no solution could 

exist any more. Thirdly, the parameter range is limited by the 
case of Se4 tending to a switch-on shock. A similar situation 
is given for the basic pattern S f 1 , T2, S3, T4 , R f 5. In the 
following we restrict ourselves to the solution involving S f 5. 

Each of the three limiting cases involves a relationship 
between the seven initial parameters. As an example we consider 
the relationship between A n , 2 and s 2 for a given set of s.,, a 2 , 

a iB,'T A n 1 and “iB 2 as shown in Figure 4. In the case A .= 8 
shown in Fig. 4a the limiting curve corresponding to the fast 
Mach number of shock S f 5 being ^ = 1 extends from s 2 = 0 to 
00 • The same is true for the limiting curve describing the tan- 
gency condition. Between these curves the solution S^I , T2 , 

Sj3, T4 , S5 is possible. Below the curve S^5 = 0 the solution 
is 3^1 , T2, Sf3, T4, R^5. The latter type of solution is limited 
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by another curve below the curve 5 = 1 expressing the tan- 
gency condition for the rarefaction wave. This curve has not 

been computed, however. In the other case A - = 4 shown in 
** n ^ 1 

Fig. 4b there must always be a reflected wave called generated 
wave by Jaggi and Wolf above 56 for a steady-state 

solution to exist in our frame of reference. Again an upper 
limit exists for the Alfven Mach number A„ , as a function of s 
The curve expressing the tangency condition extends continu- 
ously beyond the point, where it meets the curve ^ = 1 , 
thereby limiting the solutions of type S^l , T2, S^3, T4, R^5. 

Inspection of Figure 3 shows , that an upper limit to A n ^ 
due to the tangency limit can only be found if the minimum of 
curve Sf5 fulfills the inequality ot v ^ > “ « 2 ' s i nce th® curve 
S^3 is stretched progressively with increasing Mach number A^ ^ 
when all other parameters are kept constant. The marginal case 
can be expressed as a condition on , a 2 / a iB ^, an<3 A n , 1 * 

In addition to the cases illustrated in Figure 4 a large 
number of other possibilities exists. It is for example also 
possible, that only above a minimum value of M f 5 can be 
equal to one . 

The disappearance of S f 5 from the flow pattern requires 
one condition on the seven initial parameters to be fulfilled. 
S^3 cannot disappear by the argument in section II , item 5 
proving the impossibility of a wave R^3. T4 can disappear, 
however. For this to happen three conditions must be fulfilled 
by the seven initial parameters leading to s^ = s^, 

“iB ,3 - “iB ,4 and A n,3 = V-T Tha cation c iVj3 = c iVj4 
can be enforced separately by the correct choice of a iv 2 * 

A trivial case is apparently given by the parameter set 
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s 1 , ot 


2' a iB , 1 ' A n,1 ' s 2 “ S 1 ' 


a . _ _ = a . „ and A ~ = A - . 
iB,2 lB , 1 n,2 n,1 


Since the disappearance of T4 is a very special case due to 
the three required conditions on the initial parameters the 
possible existence of nontrivial solutions will not be 
considered. 


In closing this section we make two more points. In 
section Il^item 4 the assumption has been made that no switch- 
on or switch-off waves occur in the general solution. The 
development of the theory in this section is consistent with 
this assumption, since switch-on shocks occur as limiting 
cases only. Switch-off waves are not necessary at all to 
construct a solution, since they are of the S g -type. 

The second point concerns the poles in the generation 
coefficients of the paper by Jaggi and Wolf 9 . These poles do 
not occur in the nonlinear analysis. They are due to the 
linear approximations in the vicinity of points Se5 and Se4. 
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Parametric Study of the Solutions of the Interaction Problem 


In this section we present a limited parametric study of 
our interaction problem. Instead of attempting a complete 
presentation of the characteristics of the physical solution, 
we consider a few cases of geoastrophysical interest. We shall 
exclude from consideration results for large s^ and S 2 and large 
and A^ since these are well known in the gasdynamic litera- 
ture . Although the magnetic field does not greatly influence 
the results in these cases , it may be very important as a 
diagnostic tool. Here we introduce a new Alfven Mach number 


A 


A cos 
n 


iB 



(14) 


which is often more convenient than A r for applications. For 
example. 


A 1 = 


iv^n S ^ 1 j /4 tt Pi 


sin a. 


ili 


We note that the value of s = 1 is typical for the solar 
wind in the vicinity of the earth. Typical values for A 1 are 
4 and 12 for interplanetary shocks and the earth’s bow shock 
near the subsolar point, respectively. These values will not 
be chosen exclusively, however. 

The first case we consider is the interaction of S^l with 
T2, where T2 is characterized by a change in magnetic field 
direction only i.e., = p^, p^ = P 2 , v^ = y 2 or s^ = s^, 

A 1 A n,1 cos a iB,1 = A 2 = A n,2° cos a iB,2 but a iB,1 9 ener ally 
diff eren t from a ^g^ 2 * This bypo of discontinuity, corresponding 
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o 

to a type (a) discontinuity in Jaggi and Wolf's treatment of 

weak discontinuities, is most commonly found in interplanetary 
2 8 

space ’ . It illustrates most clearly the influence of the mag- 
netic field. For the typical interplanetary values of s^ = s 2 = 1 
and A.j = A 2 = 4 and a 2 = 45° Figures 5a, b illustrate various 
interesting quantities as a function of a iR ^ and a iB 2 « Figure 
5a shows Aw = W 3 ^ - e 2 -j i.e., the change of the angle between 
the B-vectors on both sides of the tangential discontinuity 


during the interaction. We have defined o>, = “'n v ’ n* 

The curves can be continued to values of a. n - > 90° by use of 
some of the symmetry properties of the solutions. Thus if a iB 2 
and a iB ^ and therefore oj^ 4 are solutions for s^ , ct 2 , a^ B 
A , s„, a._ A. „ given with a „ < 90°, say, the trans- 

n , I Z Id f z n ; Z Id / z 

formation a 1Bi2 - 180 ° - c 1B(2 yields a iB>3 - 180° - a iB>3 , 
“iB,4 * “iB,4 and therefore «> 3/4 - 180° - a iBj3 - a iB<4 and 

finally Au + 2 Am ^ a iB 2 = 90 °) “ Aa) ' i.e., 
the curves Au ( a iB 2 ^ can ke continued to <* iB 2 > 90 ° by 
reflection with respect to their points of intersection with 
the line a. n 0 = 90°. Only in the case a._. 1 = 0° an 5^5 shock 

Xd.^ Z Xd f \ X 

occurs all the time. Increasing A^ = A 2 with s^ = s 2 and a 2 
kept constant leads to a further limited increase of Aid. 


If the magnetic field could be neglected in the dynamics 


of the problem, we would expect = 0 in the case considered 
above with s^ = s 2 and A^ = A 2 . That the magnetic field influ- 
ences the dynamics is shown in Figure 5b, where angles of 
up to 4° are shown. The curves can be reflected at the line 
a iB 2 = 90° to obtain values for a iB 2 > 90°. Increasing the 
Mach numbers A^ = A 2 leads to lower values of a^, since the 
pressure forces become more and more important compared with 
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the magnetic forces. The shocks S^5 are very weak in all these 
cases. The largest value of 5 = 1.04 is attained for 
a iB 1 = °° and a iB 2 = 9°°* An interesting feature of these 
solutions is, that whereas T2 shows a change in direction of B 
(given by u> 2 -j ) only, the modified tangential discontinuity can 
display additional changes e.g. in | B j . In contrast to |B^ | = 

and a iB 2 = 90°, and a minimum density ratio of 0.9 across T4 
for the same parameters compared with a value of one at T2 . 
Increasing A^ = A 2 to a value of 8 leads to a maximum |b 3 |/|B^| 
of 1.24 at the same a . _ . and ct._ 

Id / I Id g 4 

The next case we consider involves density changes across 
T2 or equivalently a velocity shear. The results for are 
shown in Figure 6, where a 3 is plotted as a function of A 2 
with s 1 = s 2 = 1 , a 2 = + 45°, a iB 1 = ct iB 2 = 0° and A 2 = 4 and 
12. It is expected that the change in direction of the shock 
normal is much more pronounced in this case. An increase in 
density across T2 by a factor of 2, corresponding to a ratio 
A 2 /A^ = /2, leads to = - 10.9° and a fast Mach number 
r = 1.1 for A. = 4 and = - 9.3° and M- c = 1.13 for 
A^ = 12. The density ratio across T4 is 1.84 and 1.68 for 
A.j = 4 and 12, respectively, compared with two across T2. On 
the other hand the density ratios across S^5 are 1.14 and 1.21 
and magnetic field magnitude ratios 1.1 and 1.14, respectively. 

The last example illustrates the influence of a change in 
s across T2 with =0°, s 1 = 1, A 1 = A 2 = 8 on the re- 

lationship Aw («^ B 2 )° Since a change in s across T2 implies 
| B 1 | ? | B 2 | A 1 = a 2 implies an additional change in density 
or velocity v^ across T2 . The results are shown in Figure 7. 


| B 2 | we have a maximum ratio of | | / | B ^ f = 1.12 at ct iB 
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We see, that an increase in -|l implied by a decrease 

in ^ 2 ^ s, \ ^ eads to progressively smaller changes in 
Aw = u>3, 4 - a) 2J . 
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VI. Conclusions 


In conclusion we have shown that for a range of initial 
conditions the interaction between a fast shock S^I and a 
tangential discontinuty T2 can lead to a steady flow pattern 
in the frame of reference in which both discontinuities are at 
rest. For ^- n 2 — max two solutions are possible, S^I , 

T2, S^3, T4 , Sf5 and S^l , T2 , S^3, T4 , F^5, which are adjacent 
in the parameter space s^ ct 2 , A n ,i' s 2' a iB,2' A n,2* 

Uniqueness of the solutions has been enforced by the re- 
quirement of continuous connection to the limiting cases of 
linear interaction, a device, which works remarkably well in 
gasdynamics. Further work is necessary in this area concerning 
the stability of the solutions obtained. Also laboratory ex- 
periments and space observations could be used to finally 
decide the question of uniqueness. 

The set of seven dimensionless parameters turns out to be 
very economical in ordering the large number of solutions. The 
similarity laws defined thereby show that a velocity shear 
parallel to the interaction line constitutes a very weak 
interaction leading to S^3 essentially identical to S^l and no 
wave S^5 or R^5. These laws also show that velocity shears 
perpendicular to the interaction line and density changes 
across T2 are equivalent. The computations for individual 
cases show that a tangential discontinuity T2 with a change 
in direction of the magnetic field only and s^ = s 2 = 1 
produces a very weak reflected or generated fast wave, whereas 
the tangential discontinuity T4 has more complicated charac- 
teristics than T2, i.e., a change in magnetic field magnitude. 
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density, temperature and velocity shear. Because of the 

deviation of the shock normal to up to 4° for s^ = s 2 = 1 , 

<*2 = 45° and = A 2 = 4 the propagation of a shock through 

an ensemble of tangential discontinuities of this type like 

in the solar wind can lead to an appreciable random walk of 

the normal direction and other shock properties . 

An initial tangential discontinuity with a density jump 

leading to a sudden impulse in the geomagnetic field ac- 

1 5 

cordingto Burlaga splits into two disturbances after in- 
teracting with the earth's bow shock which due to their 
different propagation velocities would have different travel 
times. Since in this case S^5 or R f 5 are not weak the 
initiating signal at the magnetopause is a broad signal and 
not a one step signal as often assumed in the past. 
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appendix 


In this appendix we shall collect the functional re- 
lationships needed in equations 12 and 13. Figure 8 shows 
the geometry of the problem. The upstream region is denoted 
by the subscript zero and the downstream region by one . The 
upstream region is then characterized by s Q , a^ B Q and A n ^ Q 
and the shock or rarefaction wave by the angle a between the 
downstream edge of the Sg or Rg-disturbance and the projection 
of the streamline on the x,y-plane. This definition of a is 
useful for both Sg and Rg. For a shock the downstream and 
upstream edges are parallel trivially. For the rarefaction 
wave the direction of the upstream weak discontinuity is 
given by the initial conditions already, whereas the direction 
of the downstream weak discontinuity can be used to charac- 
terize its strength. 

In the following we restrict ourselves to the shock case. 

1 3 

We use the formulas derived by Bazer and Ericson and also 
given by Jeffrey and Taniuti^ with the evolutionarity con- 
ditions imposed by them. Note, that in our case n is pointing 
in the direction of propagation . The angle between n and B q 
is called 0 Q and without loss of generality we can write 

0 < 0 < 180° or sin 8 > 0. We obtain 

— o o — 


cos 6^ 
o 


sin a cosa 


iB ,o 


(A1) 


A shock is then possible for the angular range 

a„ < a < 180° - a M (A2) 

— — Mg 
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where cc^ is given by 


sin a 


1 


M 


£ l cos a iB ,o I s n,o 


— <‘ s o + (s o + ’> A2 n,o> V2 


(A3) 


with O < a,, < 90° and A > max ( 1 , fs ~ ) 

— n,o — ' o 


Introducing the quantity 


h = 


B. 1 - B, 
t , I t , o 

|B I 
1 — o 


(A4) 


where is the magnitude of the component of B in the shock- 
plane, we obtain the following cubic equation for h 


h 3 + h 2 \ sin e (A 2 + 11) 

3 o n,o 


+ h [a 4 4 cos 2 e + a 2 (- + 2s + 14 sin 2 e )+ f - 2s 1 

L n,o 3 o n,o 3 o 3 o 3 oj 


+ 2sin0^ I -A 4 cos 2 0 + A 2 (1 + s )-s 1 = 
o L n,o o n,o o o J 


0 


(A5) 


where we have used the ratio of specific heats X = 5/3. The 
equation has one root fulfilling the requirement 0 _< ri^ 3 
necessary for a fast shock. The relative density change is 
given by 



A ~ 1 

n ,0 


h + A sine 
n,o o 


(A6) 


The relative pressure change is given by 
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- P 1 , 5 , K! , . , n f- h sin6 o 

Y £ = p“ ' 1 = 3?~ l “ 5 h + h <- — > } 

o h - n.e sine 

f o 


We are interested in the quantities s 1( a._, A Y. ,, <5,., 

i its,! n, i t ,r r 

Using |B |/|Bl= (1 + 2 h sin 6 Q + we obtain 


- 1 + Y f 

g — — .. . . - ^ 

* s 1 + 2h sine + h 2 ° 
o 


p t,1 

_ 2 
lOn* cos 0 
t o 

^t ,o 

h-ri£ sin e Q 


o o 


A n,1 = A n,o (1 + V 


- 1/2 


(AIO) 


The quantities 6 f and a ±B 1 are a little more difficult to 
obtain: 


sin 8 q + h 

sin 8 . n 1 — sin Ct.-. _ n 4 /o t 

6,1 B ' sin 0 o (1+2h sin 0 o +h 2 ) 1/2 


(All) 


where a... is in the same quadrant as a.„ from 0° to 360° 
Id , I - 1 lB f O 

and 

2 2 2 2 ” 1/2 
sin 6 f = h sinacosa (h cos a + sin 0 +2hsin6 cos a) , (A12) 


where |(S f | < 90°. 
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Figure Captions 


Figure 1 Magnetogasdynamic steady flow pattern as a result 

of the interaction of the fast shock S^I and the 
tangential discontinuity T2 . 

Figure 2 (a) Illustration of angular ranges of "outgoing" and 
"incoming" shocks for given upstream conditions, 

(b) Illustration of "outgoing" and "incoming" fast 
rarefaction waves. 

Figure 3 Graphical representation of the uniqueness problem 
in the mathematical solution of the interaction. 
Both points Se4 as well as Se ' 4 fulfill equations 
8 and 9, i.e., P t3 = P tj4 and = a V/ „. 

Figure 4 (a) Limiting curves due to tangency condition and 

M f 5 = 1 for the steady solution S^l , T2 , £p, T4, 
S^5. For the case of R f 5 instead of S^5 only the 
limit of an infinitely weak R f 5 is shown i.e., 
the curve labeled ^ = 1 . Both limiting curves 
are completely separated in this case. 

(b) A more complicated case of limiting curves, in 
which tangency and 5 = 1 are possible at the 
same time. 

Figure 5 (a) Change Aw of the angle between magnetic field 
vectors on both sides of the tangential dis- 
continuity during the interaction as a function 
of ^ and a^ B 2 with velocity, pressure and 
magnetic field magnitude as well as density 
continuous across T2 at the beginning. 

(b) Deflection of the shock normal during the 
interaction. 
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Figure 6 


Figure 7 


Figure 8 


Deflection a 2 of the shock normal and fast mag- 
netoacoustic Mach number of the reflected wave 
Sf5 for tangential discontinuities with density 
change and/or velocity shear for = 4 and 1 2 . 

Note the large deflections compared with Figure 5a. 
A a as a function of a iB 2 for various values of s 2 
expressing magnetic field increases and plasma 

pressure decreases. T2 is of a mixed type, since 

2 


| B | changes as well as | v^ | 

Geometry of a stationary oblique shock. 
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Figure 1 
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Figure 2a 
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